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Abstract 

We study r] and r]' mesons and their mixing angle in a mixed action approach 
with so-called Osterwalder-Seiler valence quarks on a Wilson twisted mass sea. The 
gauge configurations have been generated by ETMC for Nf = 2 + 1 + 1 dynamical 
quark flavours and for three values of the lattice spacing. The main results are that 
differences in between the mixed action and the unitary approach vanish towards the 
continuum limit with the expected rate of 0{a?‘). The individual size of the lattice 
artifacts depends, however, strongly on the observable used to match unitary and 
valence actions. Moreover, we show that for the rj mass valence and valence plus 
sea quark mass dependence differ significantly. Hence, in this case a re-tuning of the 
simulation parameters in the valence sector only is not sufficient to compensate for 
mismatches in the original quark masses. 
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1 Introduction 


Mixed action approaches, where valence and sea fermion actions are chosen 
differently, are used frequently in lattice QCD. They possess a number of 
important advantages compared to the so-called unitary case, where valence 
and sea quark actions are identical. In particular, it might be possible to use 
a valence action obeying more symmetries than the sea action in cases where 
the valence action cannot be used in the sea for theoretical reasons or because 
of too high computational costs [1]. Prominent examples are overlap [2,3,4] 
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or domain wall [5] valence quarks on a Wilson-like or staggered sea. Concrete 
examples can be found for instance in Refs. [6,7,8]. 

When working in a mixed action approach, valence and sea actions need to be 
matched appropriately, for instance by tuning the valence quark masses such 
that a choice of meson masses agrees between unitary and mixed approaches. 
People even go one step further and try to correct for small mismatches in 
bare parameters used in the sea by using a partially quenched mixed action ap¬ 
proach. The most extreme example for this approach is to use valence strange 
and charm on gauge conhgurations with only Nf = 2 light dynamical quark 
flavours. In this case sea and valence actions are not matched, but the valence 
parameters are tuned such as to reproduce a choice of physical observables. 

This ansatz has the big advantage that the gauge conhgurations do not need 
to be re-generated. However, while apparently quite successful, it is question¬ 
able whether this procedure works for observables with a strong sea quark 
dependence. Due to OZI suppression there are not many examples of such 
observables. But their very existence makes a clear distinction between QCD 
and the naive quark model. 

Of course, a mixed action approach has also disadvantages, most prominently 
the breaking of unitarity, which might for instance drive certain correlators 
negative [9,10]. Also, it is not clear a priori how big lattice artifacts one en¬ 
counters in mixed formulations. 

In this paper we will present results on a particular mixed action approach 
with so-called Osterwalder-Seiler [11] valence quarks on an Nf = 2 + 1 + 
1 havour Wilson twisted mass sea [12]. This particular action combination 
has the advantage that exact valence quark havour symmetry is preserved. 
Moreover, the respective zero modes of sea and valence quarks coincide in the 
chiral limit. However, 0(0^) violations of havour (and parity) stemming from 
the sea quarks are still rehected in the magnitude of lattice artifacts on various 
physical observables. 

As physical example we study the r] and r]' system. The large mass splitting ob¬ 
served among M^o M^i is thought to be due to the Ua{X) anomaly, a 

relation established via the Witten-Veneziano formula [13,14,15]. This lets one 
expect a signihcant dependence on the sea quark degrees of freedom. Speak¬ 
ing more technically, the corresponding correlation functions obtain signih¬ 
cant contributions from fermionic disconnected diagrams and are, therefore, 
uniquely sensitive to diherences between valence and sea formulations. Note 
that this was also discussed in the context of the validity of the fourth root 
trick in staggered simulations, see Refs. [16,17] and references therein. 

After matching valence and sea actions, we compare observables extracted 
from unitary and valence operators. The unitary observables have been com- 
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puted in Refs. [18,19]. We study the continuum limit with different matching 
conditions and hnd remarkably good agreement to the unitary case. How¬ 
ever, when comparing the valence with valence plus sea strange quark mass 
dependence of we hnd signihcant differences. 

These hndings are important for future lattice QCD investigations: there 
are many phenomenologically interesting quantities involving havour singlet 
pseudo-scalar mesons, for instance form factors of B or Dg decays to r|^v. 
And maybe most prominently, there are anomaly related form factors of 
r] —)■ 77 , which can be used to estimate the light-by-light contribution to 
the hadronic part in the anomalous magnetic moment of the muon, in which 
we currently observe a deviation between theory and experiment at the few 
a level [20,21,22,23]. The usage of Wilson twisted mass fermions described in 
this paper has signihcant advantages compared to other lattice actions due to 
a powerful variance reduction. And the possibility to use a mixed action will 
further ease those computations. 

More generally, the hndings here show that with a mixed action approach 
one can deal with fermionic disconnected diagrams, provided one applies an 
appropriate matching procedure. These disconnected diagrams become more 
and more important as they need to be treated appropriately for instance 
in investigations of hadron-hadron interactions. Since we show here that a 
mixed action approach works in the case of rj, rj' mesons, where the fermionic 
disconnected diagrams contribute signihcantly, we are conhdent that the same 
approach can be used for other physical observables. First accounts of this 
work can be found in Ref. [24]. Other studies of rj and rj' mesons from lattice 
QCD can be found in Refs. [17,25,26,27,28,29]. 


2 Lattice actions 


The results we will present here are obtained by evaluating correlation func¬ 
tions on gauge conhgurations provided by the European Twisted Mass Col¬ 
laboration (ETMC) [30]. We use the ensembles specified in table 1 adopting 
the notation from Ref. [30]. More details can be found in this reference. 

The sea quark formulation is the Wilson twisted mass formulation with Nf = 
2 -|- 1 -|- 1 dynamical quark flavours. The Dirac operator for the light quark 
doublet reads [12] 

= Dw + mo + , ( 1 ) 

where Dw denotes the standard Wilson Dirac operator and /i^ the bare light 
twisted mass parameter, and in general = 1,2,3 represent the Pauli 
matrices acting in flavour space. acts on a spinor xe. = d)'^ and, hence, 

the u (d) quark has twisted mass 
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ensemble 

P 

ape 

apa 

aps 

V 

A^conf 

Nb 

A40.24 

1.90 

0.0040 

0.150 

0.190 

24^ X 48 

1117 

5 

A60.24 

1.90 

0.0060 

0.150 

0.190 

243 X 48 

1249 

5 

A80.24 

1.90 

0.0080 

0.150 

0.190 

243 X 48 

2441 

10 

A100.24 

1.90 

0.0100 

0.150 

0.190 

243 X 48 

968 

5 

A80.24S 

1.90 

0.0080 

0.150 

0.197 

24^ X 48 

2420 

10 

A100.24S 

1.90 

0.0100 

0.150 

0.197 

243 X 48 

1196 

5 

H55.32 

1.95 

0.0055 

0.135 

0.170 

32^ X 64 

4450 

5 

D45.32SC 

2.10 

0.0045 

0.0937 

0.1077 

32^ X 64 

2220 

10 

Table 1 

The gauge 

ensembles 

used in 

this study. 

For the labelling of the ensembles we 


adopted the notation in Ref. [30]. In addition to the relevant input parameters we 
give the lattice volume, the number of evaluated configurations A^conf and the block 
length Nh used for bootstrapping. Nf, was chosen such that blocks cover at least 20 
HMC trajectories of length one. 

For the heavy unitary doublet of c and s quarks [31] the Dirac operator is 
given by 

Dh = Dw + mo + • (2) 

The bare Wilson quark mass mo has been tuned to its critical value [32,30]. 
This guarantees automatic order O (a) improvement [33], which is one of the 
main advantages of the Wilson twisted mass formulation of lattice QCD. 

f] and f]' masses have been computed in this framework in Refs. [18,19,34] on 
the same set of gauge conhgurations used here (and more) - we will refer to 
this framework as the unitary approach. However, in order to account for - 
and possibly benefit from - correlations we have re-evaluated the unitary rj 
and rj' masses on exactly the same gauge conhgurations as used in the present 
study. 

The splitting term in the heavy doublet (2) introduces havour mixing between 
strange and charm quarks which needs to be accounted for in the analysis. 
However, this complication can be avoided by using a mixed action approach 
for the valence strange and charm quarks. Formally, we introduce so-called 
Osterwalder-Seiler (OS) twisted valence strange and charm quarks [11,35]. 
The Dirac operator for a single valence quark havour q reads 

Dg = Dw + mo + iUq'j^ ■ (3) 

Adapting the ideas of Ref. [11] to the 77 , 77 ' system, we introduce two strange 
and two charm quark havours, s, s' and c, c', respectively. Flavours s and s' 
will have quark mass with equal modulus, but opposite sign: fj,s = \fis\ = —Rs', 
and the same for c and c'. Formally, the lattice action is extended to include 
a fermionic action corresponding to the Dirac operators (3) for all valence 
strange and charm quark havours, accompanied by a ghost action to exactly 
cancel the contributions of the additional valence quarks to the fermionic 
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ensemble 

aMgc 

aMri^ 

aM 0 

'•ronn 

A40.24 

A60.24 

A80.24 

A100.24 

0.25884(43) 

0.26695(52) 

0.27706(61) 

0.28807(34) 

0.30708(60) 

0.31010(65) 

0.31406(46) 

0.31575(45) 

0.2375(25) 

0.2544(26) 

0.2659(25) 

0.2883(14) 

A80.24S 

A100.24S 

0.25503(33) 

0.26490(74) 

0.27168(49) 

0.27455(73) 

0.2649(16) 

0.2841(16) 

R55.32 

T)45.32sc 

0.22799(34) 

0.17570(84) 

0.26087(33) 

0.21126(34) 

0.2177(10) 

0.1494(15) 


Table 2 

Values of the unitary kaon and rjs masses Mk and in lattice units, which 
have been used to match mixed and unitary actions. In addition, we give the mass 
values of the connected neutral pion which becomes identical to rjs for mass 

degenerate light and strange quarks. This SU(3) symmetric situation is realised 
approximately for the A80.24s and yll00.24s ensembles. The kaon mass data shown 
in this table has first been published in Ref. [18]. 

determinant. For details we refer to Ref. [11]. In this reference it was also 
shown that automatic (!l(a)-improvement stays valid in this framework and 
unitarity is restored in the continuum limit. In particular, flavours s and s' (c 
and c') become identical. 

It is important to notice that at finite lattice spacing values correlation func¬ 
tions involving s and s' {fig = —f^s') differ by lattice artifacts. For instance, 
the masses extracted from the correlation function of the operator 

= ipsi'ybi’d (4) 

where the helds ipq, tjjq denote single quark helds in the so-called physical basis, 
differ from the one extracted from the operator 

(5) 

by O(a^) (we denote it with in remedy of the neutral pion in the light 
sector). Only in the continuum limit these two masses will agree again. 

Valence and unitary actions need to be matched appropriately. As shown in 
Ref. [11], in our case the matching can be performed in principle using the 
relation 

Rc/s Rct Zp jZg . (6) 

However, for the strange quark mass uncertainties in Zp/Zg are magnihed 
in aps, and thus we decided not to rely on Eq. 6. Instead, meson masses are 
used: in previous studies it was found that matching kaon masses determined 
from the operator (4) to the unitary kaon masses is best in the sense that the 
residual lattice artifacts in the results computed in a mixed action approach 
are small [36]. We will call this procedure kaon matching. 
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For details on how to compute the kaon mass in the unitary case we refer to 
Ref. [37]. We note in passing that there is no kaon mass splitting introduced by 
the twisted mass formalism in the unitary case for the choice of a degenerate 
light quark doublet |/r„| = |/id| = [32]. 

As a second matching observable for the strange quark mass we use the mass 
of the so-called r]s meson The r]s is an artihcial meson corresponding to 
the following interpolating operator 

C’S® = , (7) 

for which we assume /i* = ^s"-, unlike the s' quark considered above, which 
had opposite sign. A beneht of this particular choice is the absence of discon¬ 
nected diagrams in the corresponding two-point function. This procedure will 
be called rjg matching. For technical details, e.g. further interpolating helds 
and correlation functions we refer to section 3. 

For both matching procedures on each gauge ensemble one tunes the value of 
a^s such that the kaon or the rjg mass agrees within errors between the mixed 
and the unitary formulation. The unitary values of the masses we matched to 
are compiled in table 2. In order to compute the matching values for a/i^ we 
performed inversions on a subset of the available conhgurations in a range of 
a/is values and interpolated the squared OS meson masses linearly in a/i^. The 
matching values for a/i^ for the two matching observables and all ensembles 
can be found in table 3. The values for and at the matching points 
are compiled in the appendix in table A.6. Note that in case of matching 
we do not reach exact agreement for all ensembles within errors when recom¬ 
puting from full statistics. These numerically small differences become 
irrelevant for the rj and rj' masses themselves due to the much larger statistical 
uncertainties introduced by the quark disconnected diagrams . 

In the following we indicate quantities determined in the OS framework with 
the superscript while quantities determined in the unitary case have no su¬ 
perscript. To distinguish the two matching procedures we use the superscripts 
^ and 

As an example for the matching procedure we show in hgure 1 

and {aM^o)"^ as a function of the bare OS strange quark mass a/i^ 
for the D45.32sc ensemble. In figure 1 the aforementioned OS kaon mass 
splitting can be observed. In the limit fig = fii this splitting corresponds to 
the difference between the charged pion mass and the connected only neutral 
pion mass. The splitting is almost independent of p.*, decreasing slightly with 
increasing /i^. 

As expected, is larger than the two kaon masses and agrees with M^o in 
the limit fig = All three squared masses show a linear dependence on fig. 
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Fig. 1. (aM^®)^, as a function of the bare OS strange quark 

mass a/is for the L)45.32sc ensemble. Horizontal lines indicate the unitary mass 
values that have been used for the matching. 


ensemble 

a/i£ 

a/if 

afiT 

afJ-c 

A40.24 

0.0040 

0.02300(25) 

0.01239(25) 

0.27700(25) 

A60.24 

0.0060 

0.02322(22) 

0.01303(22) 

0.27678(22) 

A80.24 

0.0080 

0.02328(20) 

0.01338(20) 

0.27672(20) 

A100.24 

0.0100 

0.02381(21) 

0.01380(22) 

0.27619(21) 

A80.24S 

0.0080 

0.01884(16) 

0.00883(21) 

0.28116(16) 

A100.24S 

0.0100 

0.01877(22) 

0.00922(23) 

0.28123(22) 

H55.32 

0.0055 

0.01858(12) 

0.01100(10) 

0.25142(12) 

T)45.32sc 

0.0045 

0.01488(30) 

0.01180(12) 

0.17252(30) 


Table 3 

Matching values of the OS valence strange quark masses n-s for kaon and 
matching. The OS valence charm quark masses have been determined using 
Eq. (8) for kaon matching only. 

The horizontal lines indicate the corresponding unitary values that have been 
used for computing a/i^ and a/i^“. 

For the charm quark mass the estimate from Eq. 6 is less affected by uncer¬ 
tainties. In order to circumvent the need for ZpjZg^ one can re-arrange Eq. 6 
to 

fXc 2/io- fXg. (8) 

Because /is -C /ic and rj^rj' do not depend on fic, we restrict ourselves to kaon 
matching for the fig value entering the charm quark mass. The corresponding 
values for /ic = /i^ extracted in this way can be found in table 3. 

All errors are computed using a blocked bootstrap procedure to account for 
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autocorrelation as well as all other statistical correlations in the data. The 
number of bootstrap samples was taken to be 1000 and the number of con- 
hgurations per block Ni, is given for every ensemble in table 1. N/, itself was 
chosen such that the length of a block corresponds to at least 20 HMC tra¬ 
jectories of length one. This value turned out sufficient to compensate for 
autocorrelation in the observables considered in this study. 


3 Pseudo-scalar flavour-singlet mesons 


In order to extract rj and rj' states we need a set of appropriate interpolating 
operators. As we are going to work in the quark flavour basis our choice is 



X 


= IIV’c*75^c(x,f) 


in the physical basis, again denoted as V'g, With Osterwalder-Seiler valence 
fermions we have to rotate the bilinears into the so-called twisted basis denoted 
as g, g, see e.g. Ref. [38], in which also the Dirac operators in the previous 
section were written. Performing this axial rotation [12,11], one obtains the 
following operators in the so-called twisted basis 

^ ^(dii(x, t) - hM(x, t )), 

C’s(t) = -^ss(x,t), 

X 

^c{t) = - ^ cc(x, t) . 

X 

From these operators we build a correlation function matrix 

C{t)qq: = {Oq{t'+ t) , g, g G {£, s, c} , (9) 

which allows us to obtain results for masses and amplitudes; cf. section 3.2. 


The corresponding correlation functions have fermionic connected and discon¬ 
nected contributions. The case for up and down quarks is like in the unitary 
approach and discussed in detail in Refs. [39,18,19]. Therefore, we concen¬ 
trate on the disconnected contributions for strange and charm quarks. The 
correlation function of Os{t), for instance, has the following contributions 


(0,(t) O;(0))r = -TV{Gf Gf} + Tl'{Gq • T1'{G“} , 


( 10 ) 



where {.)f denotes the average of fermions only and 


G:y = {D:^){x,y) ( 11 ) 

denotes the strange OS propagator. The hrst term in Eq. 10 is the connected 
contribution and the second the disconnected one. Note that mixed flavour 
correlation functions have only disconnected contributions by dehnition. The 
ground state mass extracted only from the connected piece on the r.h.s. of 
Eq. 10 is the mass of the artihcial r]s meson, which is employed for rjs matching. 

We evaluate the connected only contribution to Eq. 10 using the one-end- 
trick [40]. In contrast to the Wilson case, TrlG^^G**^} is in general complex 
valued. However, the imaginary part of the corresponding trace is a pure lattice 
artifact. This can be shown by considering a suitable combination of connected 
correlation functions involving OS quarks s and s' 

- {s -t s'})f = (ss(a;)ss(0) - s's'(a;)s's'(0))i;’ 

= - Tr{G°"Gf} + TiiG^JGf} 

= - Tr{G°"75(G°f)t75} + Tr{G°f75(G°")H5} 

= - Tr{G°"75(G°f)t75} + TT{Gf^,{Gfy-f,y 
= -2^ImTr{G°"75(G°f)t75} 

= -2i ImTr{G°^Gf}, 

where we have used the relation Dg = 75 / 1^,75 together with the cyclic prop¬ 
erty of the trace. Since the l.h.s. of the above relation vanishes in the continuum 
limit, we will drop the imaginary part in our calculations. 

For the disconnected contribution to Eq. 10 we need to estimate Tr{G**} on 
every gauge conhguration and all t-values. Tr{G**} is again in general complex 
valued. And again, one can show that the real part is a pure lattice artifact. 
Similar to the case of the connected contribution this can be inferred from the 
following equality 

{-^shb^six) + i’s'i''f5^s'{x))F = {ss{x) + s's'{x))f 

= -Tr{Gf } - Tr{G^f } 

= -Tr{Gr} - TrlGf}^ 

= - 2 ReTr{Gf}, 

which is zero in the continuum limit. Therefore, we will also drop the real part 
of the disconnected loops in the calculation. Similarly one can show that 

{ss{x) — s's'{x))f = —ImTrjGg"^} . ( 12 ) 

We remark that all of the above results hold for any further valence quark as 
well (e.g. the charm quark). 
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The full strange correlation function after subtraction of lattice artifacts is 
then given as 

(a(t) OKO))f = -ReTr{G°'Gf} - ImTrlG"} ■ ImTr{Gf} , (13) 

and analogously for the charm. Cross flavour terms involve only disconnected 
diagrams and are for instance given as 

(a(t) OI{0))f = -ImTrlG**} ■ ImTrlG^} . (14) 


3.1 Variance Reduction 


The relation (12) enables us to use a very powerful variance reduction method 
developed originally for the disconnected contributions of the light doublet [39] 
also for strange and charm flavours (see also Ref. [41]). It is based on the 
identity (recall = —/i^') 

= -2i^isD-} 75 D-^. 

Therefore, using Eq. 12 we can estimate 

ImlV{Gf} = -li.-Irla? Gf} , (15) 

and correspondingly for the charm quark. Following Ref. [39], we apply this 
variance reduction method also to the light doublet. 


3.2 Matrix of Correlation Functions 


By applying these results, we compute the matrix of Euclidean correlation 
functions in Eq. (9) and solve the generalised eigenvalue problem (GEVP) 
[42,43,44] 

Cit) v^^\t,to) = A(")(t,to) Cito) n(")(t,to), (16) 

for determining the meson masses M^i (and possibly Mr^f) from the prin¬ 
cipal correlators \^"'\t,to), n G ri,ri'. The effective masses are then computed 
by numerically solving 

(t, to) exp“^eff^* -)- exp“^iff^ 

\C)(t-\- l,to) 4- 

for aMgg^. The matrix C is enlarged to a 6 x 6 matrix by using in addition 
fuzzed [45,40] operators. 
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At this point we recall that in the unitary case also the mass of the rjs is 
not obtained from a single correlation function but rather from the ground 
state of a correlation function matrix involving connected correlation func¬ 
tions for strange and charm quarks. This minor complication arises due to 
the violation of flavour symmetry in the Wilson twisted mass formulation 
and the fact that the action can no longer be chosen flavour-diagonal for a 
non-degenerate doublet. Therefore, one has to consider off-diagonal connected 
correlation functions in addition to the ones consisting only of strange and 
charm quarks. However, the off-diagonal connected contributions are a pure 
lattice artifact and in the continuum limit the expected behaviour is restored, 
i.e. strange and charm sector decouple regarding the connected pieces. 

Apart from meson masses also matrix elements can be extracted from the 
GEVP, which are needed to obtain 77 and r]' mixing angles. We dehne the 
mixing angles (j)£, 0 * in the quark flavour basis using the pseudoscalar density 
matrix elements Aq^n = ( 0 |Oq|n) with n G {rj, t]'} and q G {/, s} as 

Ae,r] ^ /qcos 0 £ -c^ sin 0 s 

Ae^n' As^n'j \ci sin (j)i cos 0 s 

see also Refs. [18,19]. From chiral perturbation theory combined with large Nc 
arguments |0f—0s|/|0£+0s| ^ 1 can be inferred according to Refs. [46,47,48,49] 
which is conhrmed by lattice QCD [19]. Therefore, we will consider only the 
average mixing angle 0 

tan ^0 = - ^^’"f^’V (18) 

Ai^nAs,ri' 



3.3 Excited State Removal 


To improve the r]' (and rj) mass determinations, we use a method first proposed 
in Ref. [50], successfully applied for the 772 (the rj' in Nf = 2 flavour QCD) in 
Ref. [39] and very recently to the Nf = 2 + 1 + 1 case in Ref. [19]. It relies on 
the assumption that disconnected contributions are signihcant only for the 77 
and 77 ' state, but negligible for higher excited states. This means, in turn, that 
only the connected contributions to C are affected by excited states. 

Since the signal-to-noise ratio of the connected contributions is much larger 
than the one of the disconnected ones, we can determine the corresponding 
ground state at large Euclidean times and subtract the excited states at small 
times. This subtracted connected and the full disconnected contributions are 
combined in C®“’?(f), which is then used in the analysis. We refer to the dis¬ 
cussion in Ref. [19] for more details. 

This procedure clearly depends on the validity of the assumption. However, it 
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can be checked in our Monte-Carlo data: if the subtracted combination of con¬ 
nected and disconnected contributions does not show excited states anymore, 
we take it as a strong hint for the validity of the assumption. This is the case 
for all our ensembles, and it was also the case in the unitary approach [19]. 

For all ensembles we observe that the mass of the rj meson is unaffected by 
this procedure within errors. Only the error estimates get signihcantly smaller. 
This is also the case for the r]' where, however, the errors are quite signihcant 
before excited state removal. 

As an example we show in hgure 2 the effective masses of the two lowest-lying 
states for the A80.24s ensemble from r]s matching. The values shown in the 
left panel are obtained from the standard method using a 6 x 6 correlation 
function matrix build from local and fuzzed operators, whereas the right panel 
shows the results for a 3 x 3 local-only correlation function matrix with excited 
states removed from its connected contributions. Our htted values are always 
obtained from a cosh-type £t to the respective principal correlators. The cor¬ 
responding £t ranges are indicated by the bands in the plots. In general, our 
choice for the £t ranges [^ 1 ,^ 2 ] and the values of to for the GEVP are given 
in table A.7. We remark that - since there is usually no clear plateau reached 
for the 7]' state from the standard method - we apply a two-state £t to the 
corresponding principal correlator in this case. In all other cases we employed 
a single state fit in the plateau region. 

Comparing the two panels in hgure 2 one observes that the r] mass plateau 
is unaffected, but starts at t/a = 2 in the case of removed excited states. For 
the f]' there is no plateau reached in the left panel before the signal is lost in 
noise, whereas in the right one a reasonable plateau is visible. The extracted 
masses still agree within errors. 


4 Results 


In order to compare the mixed case with the unitary case we match the two 
actions as detailed in the previous sections using either the kaon or the r]s 
mass. Next we compute OS meson masses at these matching points. As an 
example we show in figure 3 the effective masses for the principal correlators 
to) of Tj and 77 ' as a function oit/a after removal of excited states from 
the connected contributions. For rj and r]' a plateau in the effective masses 
is visible from early t/a on. The corresponding result of an exponential fit is 
indicated by the horizontal lines. The £t range corresponds to the extension 
of the lines in t/a. All OS meson masses are compiled in the appendix in 
tables A.l and A. 2 . 
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Fig. 2. Effective masses and for the ^80.24s ensemble using matching 
for (a) a 6 X 6 correlation matrix including local and fuzzed operators and (b) a 3 x 3 
local-only correlation matrix with connected excited states subtracted. The fitted 
values are shown as lines with error band. The corresponding fit range is indicated 
by the length of the lines. For further details see text and table A.7. 



(a) (b) 


Fig. 3. We show the effective masses and from a 3 x 3 correlation 

matrix after subtraction of excited states as described in the text for ensemble 
AGO.24. Panel (a) is for kaon matching and (b) for matching. 

The ?7c decouples from ly and ly' and the corresponding signal is lost in noise 
very early in t/a. Hence, we will not discuss it further here and due to the 
decoupling we will also not discuss the charm quark mass dependence of op¬ 
erators in the following. 

It turns out that the choice of the matching variable makes a signihcant dif¬ 
ference for the extracted value of M^. Moreover, we always hnd < aM^. 
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tja tja 

(a) (b) 

Fig. 4. We show the effective masses aM^ and of the two lowest-lying states on 
the B55.32 ensemble for (a) kaon matching and (b) the unitary case using subtrac¬ 
tion of excited states as described in the text. 

On the other hand, the value of is unaffected within statistical errors. We 
hnd this consistently for all the ensembles investigated; cf. tables A.l and A.2. 

In addition one observes < (jB’’ by approximately 15°; cf. table A.3. This 
results from a change in the overlap of mass and flavour eigenstates, leading 
to an increased light quark contribution to the rj for kaon matching. Conse¬ 
quently, the light quark contribution to the r]' is reduced, while the respective 
strange quark contributions behave in the opposite way. Since most of the 
noise is introduced by the light quark disconnected diagrams in our calcula¬ 
tion, this explains why in general exhibits a smaller statistical error than 
whereas the error for is larger than the one for . 

However, there is a tendency that kaon matching leads to worse plateaus 
than rjs matching. A particularly extreme case of this behaviour is shown in 
hgure 4 for the H55.32 ensemble. In the left panel the effective masses for the 
two lowest lying states from kaon matching are plotted. Clearly there is no 
plateau visible for the hrst excited state. For comparison and to guide the eye 
we show the situation in the unitary setup in the right panel, calculated on 
the same set of conhgurations. This is the only ensemble for which we cannot 
identify a plateau for the rj' safely. Therefore, we will not quote a value for the 
7]' mass for 555.32 and kaon matching. 

Although the observed behaviour on 555.32 can still be interpreted as a statis¬ 
tical fluctuation, it might - in principle - also be caused by unitarity violation. 
However, it is neither possible to verify nor exclude the latter from our present 
data. While earlier studies [6,51] observed a sign flip in the scalar correlator 
signalling unitarity violation at least for a certain regime of valence quark 
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masses, a similar argument cannot easily be extended to our case. The reason 
is that only the strange qnark is treated in a mixed action approach while 
the light quarks are unitary. When looking at the scalar correlator made from 
strange qnarks only, we do not observe it to become negative on any of our 
ensembles. 

Another observation regarding the two matching methods concerns the be¬ 
haviour of the ground states in the correlation functions used to build the 
full correlation fnnction matrix. One expects all correlators with the same 
quantum numbers to asymptotically approach the same ground state mass. 
We observe this for the nnitary case, where the rj mass can be extracted from 
all correlators in the matrix C (diagonal and off-diagonal) at large Euclidean 
times. In the OS case we observe a similar behaviour for rjs matching, bnt for 
kaon matching e.g. the strange-strange correlator alone does often not repro- 
dnce the rj mass from the light-light correlator. This might signal nnitarity 
violations for the kaon matching procednre on the one hand, and can explain 
the worse plateaus for this particular choice of the matching observable on the 
other hand. 

Finally, for we observe the error to be reduced approximately by a factor 
of two with respect to the unitary resnlt for all ensembles (cf. table A.l). 
We attribnte this to the fact that we can use the variance rednction method 
discussed in section 3.1 for the OS strange quark, which is not possible for the 
unitary strange qnark. However, the errors of and do not show such 
an error reduction (cf. table A.2 and A.3), presnmably because the strange 
quark contributes little to these observables. 


4-1 Light Quark Mass Dependence 


The first goal of this paper is to compare unitary to mixed action approaches 
and study the continuum limit of the corresponding differences. For this pur¬ 
pose we will study differences of qnantities of the form AO = 

Our ensembles at different valnes of the lattice spacing are not at exactly 
identical light and strange sea qnark masses. Therefore, we have to under¬ 
stand whether we can nevertheless stndy the continunm limit. 

Theoretically, the answer to this question is yes: both in the nnitary and in 
the mixed action approach we may write 

Qiat _ + C>(a2A|cD) 

and henceforth the aforementioned difference AO is 0{o?‘A^QYi)- ^ quark mass 
dependence is expected to be negligible because for the quark mass difference 
djj, Aqcd holds. Since AO is always computed on identical gange conhg- 
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Fig. 5. (a) roAM^ as a function of (roMps)^ for the four A ensembles A40, AGO, 
A80(s) and AlOO(s). Open symbols correspond to the s-ensembles. (b) like (a), but 
for roAMfji. The kaon matching and s-ensemble data have been displaced horizon¬ 
tally for better legibility. 


/3 

1.90 

1.95 

2.10 

ro/a 

5.31(8) 

5.77(6) 

7.60(8) 


Table 4 

The chirally extrapolated values for ro/a at each value of {3 corresponding to the 
three different lattice spacings [8]. 

urations there is no physical quark mass dependence that needs to be taken 
into account, because it cancels in the difference. 


Despite this theoretical argument, let us also investigate this point numerically. 
We hrst study the light quark mass dependence of the difference between 
unitary and OS values of Mn and M^/. For this purpose we focus on the A- 
ensembles A40, A60, A80 and AlOO, where we have different light quark mass 
values available. We denote 

AMx = , X = r],r]' 

the difference between unitary and OS meson masses. Analogously we dehne 
the angle difference A(j). r^AM^ is shown for the A-ensembles in the left panel 
of hgure 5 as a function of (roMps)^. The chirally extrapolated values of rg/a 
used in this study have been determined in [8] and are listed in table 4. Filled 
circles represent the rjs matching results, hlled boxes the corresponding kaon 
matching results. The differences are computed using exactly the same con- 
hgurations leading to reduced statistical errors due to the correlation between 
unitary and OS data. 

For T]g matching r^AM^ is for all four investigated ensembles compatible with 
zero within one sigma, while for kaon matching the difference is always positive 
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and not compatible with zero. For both matching procedures, but in particular 
for rjs matching, the dependence on (roMps)^ is not signihcant within our 
statistical uncertainties. 

In the right panel of hgure 5 we show as a function of (roMps)^. 

Despite the larger uncertainties, the differences are compatible with zero for all 
ensembles and both matching procedures. There is a slight trend for differences 
with larger modulus for kaon matching. Like for the f] the light quark mass 
dependence is not signihcant. 

The angle difference A0 shows a very similar behaviour to AM,^, see tables A.4 
and A.5. For rjg matching the difference is compatible with zero, while for kaon 
matching a value of about —15° is observed. Also here the light quark mass 
dependence is not signihcant. Besides, we hnd that the diherence between (j)£ 
and (ps is compatible with zero for both matching methods and compatible to 
the one found in the unitary setup [19], again conhrming the smallness of OZI 
suppressed corrections. 

In order to check the strange quark mass dependence of the diherences AM 
and A(j) we make use of the A80, A80s and AlOO, AlOOs ensembles. The s- 
ensembles diher from their non-s counterparts only by a diherent bare strange 
quark mass value. The corresponding values for the diherences dehned before 
are also displayed in hgure 5 with open symbols. For rjs matching the diher¬ 
ences show no dependence on the strange quark mass, whereas this cannot 
be concluded completely for kaon matching. In particular, we see for AlOOs 
deviations for kaon matching, but statistical errors can still account for the 
deviation. 


4.2 Continuum Limit 


Next, we study the dependence on the lattice spacing. For this purpose we 
use the ensembles A60.24, B55.32 and Z145.32sc, which have approximately 
the same physical value of the pion mass, i.e. roMpI^-^"^ = 0.917(14)stat, 
roMpg^^'^^ = 0.888(09)gtat and = 0.911(ll)stat, where we included 

the statistical error from the respective, chirally extrapolated values of rg/a. 
As discussed in the previous section, we do not expect the residual diherences 
in the light and strange quark masses to have any ehect on this study. 

The diherence between OS and unitary results AM for both matching proce¬ 
dures is shown for r] and r]' in hgure 6 as a function of {a/r^Y ia fhe left and 
right panel, respectively. The lines represent linear hts in {a/r^Y to our data, 
and the corresponding continuum extrapolated values are shown with open 
symbols. 
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(a) (b) 

Fig. 6. Continuum extrapolation of (a) and (b) rgAM^/ as a function of 

(ro/a)^. We show the continuum extrapolated results from a linear extrapolation of 
the three ensembles (D45.32sc, i?55.32, ^60.24) in (a/ro)^ as open symbols. The 
continuum results are displaced horizontally for legibility. 


is shown in the left panel of fignre 6. For both matching procednres 
we observe a linear dependence in (a/ro)^. A corresponding continnnm ex¬ 
trapolation in (a/ro)^ leads to the expected vanishing of this difference at 
a = 0 within errors. Kaon matching clearly exhibits larger a} artifacts, while 
Tig matching gives tqAM^ compatible with zero for each value of the lattice 
spacing separately. 

In the right panel of hgure 6 we show AM^/, again for both matching pro¬ 
cedures. We remark that for kaon matching it is not possible to perform a 
£t from our present data, because of the missing mass value on the 555.32 
ensemble which is due to a bad plateau, as discussed above. In this case, sta¬ 
tistical errors are signihcantly larger. However, within their larger errors the 
difference for the two matching procedures seems compatible and the differ¬ 
ence vanishes in the continuum limit for rjg matching, as indicated by the fitted 
line in the plot. In contrast to the rj mass, it cannot be concluded that lattice 
artifacts for kaon matching are signihcantly larger than for rig matching. 


In the left panel of hgure 7 Acf) is shown as a function of {a/roY, again for 
the ensembles A60.24, 555.32 and 545.32sc. Like for AM^ we observe also 
for A0 larger diherences for kaon matching compared to rjg matching. For rjg 
matching the diherence is in fact compatible with zero for all three ensembles 
separately. For both matching procedures the continuum extrapolated values 
are compatible with zero. 
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(a) (b) 


Fig. 7. (a) like figure 6, but for the mixing angle difference Ac/). In (b) we show 
at 2 GeV in the MS scheme as a function of (a/ro)'^ for the two methods Ml 
and M2 to estimate Zp presented in Ref. [8]. 

Finally, we show in the right panel of figure 7 the quark mass difference 

= ( 19 ) 

Zp 

as a function of (a/ro)^ for the three ensembles ^460.24, B55.32 and D4:5.32sc 
at 2 GeV in the MS scheme. The renormalisation constant Zp has been taken 
from Ref. [8]. The two colours correspond to the methods Ml and M2 for 
estimating Zp. We refer to Ref. [8] for the details. In the continuum limit it 
is expected that the two matching conditions agree and the difference should 
vanish like a^. This is what is confirmed by figure 7 (b). Also, the two methods 
Ml and M2 give compatible results in the continuum limit, as expected. 


4-3 Dependence on Sea and Valence Strange Quark Mass 


Next we study the dependence of the r] (and in principle also the p') meson 
mass on the valence and sea quark mass values. As said in the introduction, 
the dependence on the valence and sea quark masses must be identical (at 
least within errors) to legitimate re-tuning in the valence quark masses only 
against sea strange quark mass mismatches. For this purpose we first define 
the dimensionless quantity 


fixed sea ensemble 

which can be computed using the two matching points we have available for 
each ensemble. For estimating from two values at each ensemble, we 





D = 

^ d(MOS)2 
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ensemble 

2^val 

7? 

A40.24 

0.76(22) 

A60.24 

0.97(17) 

A80.24 

1.17(17) 

AlOO.24 

0.92(18) 

A80.24S 

1.07(11) 

AlOO.24s 

1.32(08) 

S55.32 

0.93(09) 

Z145.32SC 

0.86(31) 


Table 5 

The valence derivative -D™* obtained in the OS case by using the mass values from 
the kaon and rjg matching points. 


ensemble 


DOS 

v 

A80/A80S 

1.54(13) 

1.37(07) 

AlOO/AlOOs 

1.34(15) 

1.67(11) 


Table 6 

We list the values for evaluated both for the unitary and the OS case using the 
two A80 and the two AlOO ensembles. This derivative includes, in contrast to D™*, 
the valence and sea strange quark mass dependence. The values for are for r]s 
matching 



(aA/OS)2 

Fig. 8. (M0S)2 as a function of for ensembles A80.24 and A80.24s. 

have to assume that depends linearly on to a good approxi¬ 

mation. That this is the case can be seen in hgure 8, where we show 
as a function of for the two ensembles A80.24 and A80.24s. The lines 

represent linear hts to our data, which describe the data well within errors. 
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We expect D™* to be mostly sensitive to the valence strange quark mass if 
computed for several valence p-s-values on the same ensemble. The results for 
are compiled in table 5 for all ensembles and r]s matching. They appear 
- independently of the lattice spacing, light and strange quark mass values - 
to be all compatible with 1. Taking the weighted average we obtain D™* = 
1.09(5), where the error is purely statistical. Concerning possible systematics 
we stress that there is no trend visible from the data, e.g. regarding a quark 
mass or lattice spacing dependence. 

Including the sea strange quark mass dependence, the corresponding derivative 
is given by 

_ d{M,f 

" d{MKf ■ 

For the unitary case we hnd = 1.45(10) and in the OS case = 1.46(6), 
using the AlOO.24, Al00.24s and A80.24, A80.24s ensembles and rjg matching, 
see table 6. While Dn and are compatible within errors, they differ sig- 
nihcantly from For kaon matching the relative statistical errors on 
turn out to be at least a factor hve larger than for rjg matching. Therefore, a 
meaningful statement regarding the compatibility of with and 
for kaon matching is not possible from our current data. 

We take this as an indication that is indeed a quantity with a signihcant 
sea strange quark mass dependence. Therefore, correcting for mismatches of 
the sea strange quark mass value in the valence sector only is not enough for 
M^. 

In principle, the difference between and that we found for a single lat¬ 
tice spacing could also be a lattice artifact. We do not think this is the case for 
two reasons: hrst, in Ref. [18] the value of was used to correct a mismatch 
in the strange quark mass tuning for all three lattice spacings available. And 
we did not observe large cut-off effects introduced by this procedure. Second, 
also is merely independent of the lattice spacing. In fact, as the Dx are 
computed from differences and the leading lattice artifacts are independent of 
the quark mass, it is expected that these quantities are not plagued by large 
cut-off effects. 

Unfortunately, the statistical uncertainty on the r]' meson masses is too large 
to allow for a meaningful investigation of Within errors this quantity is 
always zero, irrespective of whether the valence or the full strange quark mass 
dependence is considered. Moreover, has a larger light than strange quark 
contribution. It would, therefore, be interesting to perform the same study for 
the valence light quark mass instead of the valence strange. 

Finally, we remark that it is in principle possible to calculate the difference 
between and D™* from chiral perturbation theory [52,53,54]. At leading 
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order the corresponding prediction is 4/3 for both derivatives, implying that 
the difference is an NLO effect. 


5 Summary and Discussion 


In this paper we have studied rj and rj' mesons in a mixed action approach 
and in comparison to the unitary results. The mixed action was so-called 
Osterwalder-Seiler fermions on a twisted mass sea with Nf = 2-I-1-I-1 dy¬ 
namical quark flavours. 

We have found that indeed the difference between mixed and unitary results 
vanishes as the continuum limit is approached. The rate is as expected of 
0{a?) [ 11 ] for all quantities and matching procedures investigated in this pa¬ 
per. 

For the rj mass we hnd a signihcant dependence of the size of the cutoff effects 
on the matching procedure. Lattice artifacts in the difference to the unitary 
result are compatible with zero when the two actions are matched using the 
rjs meson, while they are of normal size when the kaon is used as a matching 
variable. The same is true for the mixing angle and the strange quark mass. 

For Mjji we do not observe a strong dependence on the matching procedure. 
This can have two reasons: first the error of M^/ is large making precise state¬ 
ments difficult. Second, the r]' receives a strong contribution from sea quarks, 
because it is mainly the singlet state. The sea quark contributions are un¬ 
affected by different choices of the valence strange quark mass. Hence, this 
hnding might reflect the physical properties of the f]' meson. 

This shows that the mixed action approach can also be applied in practice 
for flavour singlet quantities and, more generally, for observables involving 
fermionic disconnected diagrams. In case of r]s matching we hnd even reduced 
statistical errors for the r] mass which might turn out to be an important ad¬ 
vantage of the mixed approach. Thus, we will use the mixed action to inves¬ 
tigate more complicated problems like ?7 —)■ 77 form factors or Ktt scattering 
for / = 1/2 in the future. 

Another important result of this paper is that for the rj meson it is not suf- 
hcient to re-tune the valence quark masses to correct for small mismatches 
in the simulation runs. Our data shows that the valence strange quark mass 
dependence of the rj differs signihcantly from the dependence on the sea plus 
valence strange quark mass. And this difference is not vanishing as the con¬ 
tinuum limit is approached. For the 7 ' we cannot make such a statement due 
to too large statistical errors, but we expect a similar result once M^/ can be 
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determined with higher accuracy. 


The latter hnding seems to contradict previous studies where the valence quark 
masses have been re-tuned to their physical values in the same mixed action 
approach for instance to determine non-singlet pseudo-scalar decay constants 
or quark masses (see for instance Refs. [8,55,56]). However, these investiga¬ 
tions were concerned with observables for which the quark mass dependence 
is expected to be mainly governed by the valence quarks. For the rj and rj' 
mesons studied here this is not the case as OZI violating contributions are 
anomalously large. However, with high enough accuracy the effect seen here 
should also show up in other physical quantities, but on the current level of 
precision it is likely to be a negligible systematic uncertainty. 
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A Data Tables 


In this appendix we have compiled all data in tables for convenience. 


ensemble 

aMrj 


aM^ 

A40.24 

A60.24 

A80.24 

A100.24 

0.2837(47) 

0.2870(28) 

0.3009(20) 

0.3074(23) 

0.2793(35) 

0.2899(16) 

0.2970(09) 

0.3067(13) 

0.312(12) 

0.3285(74) 

0.3410(67) 

0.3412(75) 

A80.24S 

A100.24S 

0.2678(13) 

0.2759(17) 

0.2690(07) 

0.2741(10) 

0.3130(42) 

0.3247(30) 

.655.32 

T)45.32sc 

0.2467(12) 

0.1890(34) 

0.2465(07) 

0.1929(39) 

0.2753(34) 

0.2032(72) 

Table A.l 




Results for aM„ 

for the unitary and the mixed action approach. For the latter we 

show the values corresponding to rjg and kaon matching. 


ensemble 

aMrji 

aM^>; 

aMp 

A40.24 

A60.24 

A80.24 

A100.24 

0.443(27) 

0.482(27) 

0.481(26) 

0.461(24) 

0.457(28) 

0.474(27) 

0.485(29) 

0.441(21) 

0.448(15) 

0.458(15) 

0.466(17) 

0.442(13) 

A80.24S 

A100.24S 

0.465(23) 

0.542(42) 

0.461(25) 

0.523(40) 

0.431(13) 

0.463(20) 

R55.32 

D45.32SC 

0.425(12) 

0.278(12) 

0.415(12) 

0.269(12) 

bad plateau 

0.271(09) 

Table A.2 




Same as table A.l, but for aM^/. 



ensemble 

4> 



A40.24 

A60.24 

A80.24 

A100.24 

47.1(2.0) 

49.2(1.4) 

49.9(1.1) 

49.6(1.3) 

47.3(2.2) 

49.8(1.6) 

51.1(1.1) 

50.3(1.3) 

30.6(4.1) 

34.2(3.7) 

38.0(3.5) 

31.9(4.1) 

A80.24S 

A100.24S 

51.3(0.8) 

53.8(0.7) 

52.6(0.5) 

55.1(0.6) 

36.2(3.2) 

44.2(2.7) 

R55.32 

Z145.32SC 

48.2(0.8) 

45.6(2.5) 

49.5(1.3) 

47.3(4.4) 

34.0(2.9) 

38.2(5.8) 


Table A.3 


Same as table A.l, but for the mixing angle (/>. 
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ensemble 

aAMrj 

aAMrj' 

A(p[Deg] 

A40.24 

+0.0281(89) 

+0.004(14) 

-16.5(2.6) 

A60.24 

+0.0415(61) 

-0.024(15) 

-15.0(2.7) 

A80.24 

+0.0401(58) 

-0.015(11) 

-11.9(2.6) 

A100.24 

+0.0338(69) 

-0.018(13) 

-17.7(3.3) 

A80.24S 

+0.0452(43) 

-0.034(12) 

-15.1(2.8) 

A100.24S 

+0.0489(31) 

-0.079(24) 

-9.5(2.7) 

555.32 

+0.0286(35) 

N/A 

-14.3(2.9) 

545.32SC 

+0.0142(44) 

-0.007(06) 

-7.4(3.7) 

Table A.4 




Values for aAM^, aAM^ 

! and A(f)[Deg\ 

as defined in the main 

text for kaon matching 


procedure. 


ensemble 

aAMri 

aAMrji 

A(t>[Deg\ 

A40.24 

-0.0044(32) 

+0.014(11) 

+0.2(1.2) 

A60.24 

+0.0029(24) 

-0.008(10) 

+0.5(1.0) 

A80.24 

-0.0039(18) 

+0.004(08) 

-1.2(0.6) 

A100.24 

-0.0007(22) 

-0.020(09) 

+0.7(1.0) 

A80.24S 

+0.0011(15) 

-0.004(08) 

+1.3(0.9) 

A100.24S 

-0.0018(20) 

-0.019(13) 

+1.3(0.8) 

555.32 

-0.0002(14) 

-0.010(17) 

+1.3(1.5) 

545.32SC 

+0.0039(20) 

-0.010(05) 

+1.8(2.4) 


Table A.5 

Same as table A.4, but for r]s matching procedure. 


ensemble 

aM^ 


'Is 

aM^; 

A40.24 

0.2583(16) 

0.2031(17) 

0.3824(13) 

0.3046(18) 

A60.24 

0.2669(09) 

0.2161(09) 

0.3838(10) 

0.3123(14) 

A80.24 

0.2770(06) 

0.2297(06) 

0.3860(08) 

0.3126(12) 

A100.24 

0.2880(08) 

0.2422(09) 

0.3923(11) 

0.3178(16) 

A80.24S 

0.2551(07) 

0.2027(07) 

0.3533(09) 

0.2720(12) 

A100.24S 

0.2649(07) 

0.2173(07) 

0.3538(12) 

0.2695(18) 

555.32 

0.2280(04) 

0.1893(04) 

0.3221(04) 

0.2640(06) 

545.32SC 

0.1758(10) 

0.1618(10) 

0.2319(04) 

0.2096(05) 


Table A.6 


OS kaon and rjs mass values for both matching procedures. 
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ensemble 




fils 



^2 



^2 

A-ensembles 

12 

22 

12 

22 

1 

3 

12 

1 

2 

5 

555.32 

15 

25 

15 

25 

1 

3 

16 

1 

2 

5 

545.32SC 

18 

30 

18 

30 

1 

3 

16 

1 

2 

5 

A80.24^^^ 





2 

7 

15 

1 

2 

8 


Table A.7 


List of fit parameters. and define the fit intervals for the ground state of 
the connected correlation function in the light (strange) sector, which is required for 
the subtraction of excited states from the full correlator, before solving the GEVP 
for the resulting correlation function matrix. The to values for the GEVP used to 
determine the r], r( states are given by whereas denote the respective 

fit ranges to the principal correlators. In the last row we give the parameters for the 
GEVP applied to the full 6x6 correlation function matrix of the A80.24 ensemble 
for the ijs matching case, as shown in the left panel of figure 2. 
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